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1 Introduction 



Vacuum decay is an old subject in field theory [1]. Coleman and Callan [2] showed that a 



quantum tunneling process from a false vacuum to a true vacuum can be realized via the 
nucleation of a true vacuum bubble in the surrounding of a false one. Coleman and De 
Luccia [3] found that gravity has a significant effect on the vacuum decay process. 

In semiclassical approximation, the decay rate per unit volume is given by an expression 
of the form 

r = Ae~ SE/h [l + 0(%)\, (1) 

where Se is the Euclidean action for the bounce: the classical solution of the equation of 
motion with appropriate boundary conditions and the prefactor A, comes from Gaussian 
functional integration over small fluctuations around the bounce, and has been discussed in 
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Ref. [3, 4, 5]. The solution of the equation of motion looks like a bubble in four dimensional 
Euclidean space with radius R and, thickness proportional to the coefficient of the symmetry 
breaking term in the potential. If there is more than one solution satisfying the boundary 
conditions, then that with the lowest Se dominates Eq. (1). 

Recently, some authors have discussed the vacuum decay in different situations such as: 
different scalar field theories [6], nonminimal coupling between the scalar field and curvature 
scalar [7], DBI action [8], and non-thin- wall limit [9], etc. The finite temperature effect on 
the false vacuum decay process has also been discussed by Linde et. al. [10], where one 
should look for the 0(3)-symmetric solution due to periodicity in the time direction with the 
inverse temperature period T _1 , instead of the 0(4)-symmetric solution at zero temperature. 
The cosmological applications of false vacuum decay process have been applied to various 
inflation cosmological models [11]. 

In general, it is not possible to find an analytical solution of the field equations for a 
finite difference of the potential between two minima. An approximation, where the field 
equations can be solved exactly, is the "thin wall approximation" (TWA) [2]. An analytical 
calculation of the nucleation rate for first order phase transitions beyond the TWA of the 
standard Ginzburg-Landau potential with asymmetric term were studied by Minister and 
Rotsch [5]. In this paper, we extend our earlier work for the theory with 3 symmetry 
breaking term [6], and we calculate the bounce and the radius of the bubble up to second 
order by expanding the bubble solution in powers of the asymmetry parameter. 

The paper is organized as follows. In Sec. 2 we present the Euclidean action and the 
equation of motion of the scalar field 0. In Sec. 3 and 4 we calculate the bounce and radius 
of the bubble as well as the action. Sec. 5 includes our conclusion and discussion. 

2 Equation of the motion and the action 

Let us consider a three-dimensional scalar field theory with a Lagrangian density 

£=^0^0-f/ s (0), (2) 

where U s ((f>) is a symmetric double- well potential having two degenerate minima at = ±1 
and has the from 
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Us(<P) = l(<P 2 -l) 2 - (3) 
To shift mutually the minima, we introduce a small asymmetry term proportional to 5: 

f/(0) = i(0 2 -l) 2 + 50 3 , (4) 

where < 5 < 2. The parameter 5 fixes the asymmetry of the potential. In particular, the 
difference between the values of the potential (f>± is 

U{4> + ) - C/(0_) = 25 + 0(5 3 ) (5) 

The shift of the minima in terms of the asymmetry parameter 5 is given by 

± = ±l-^±^ 2 + O(<5 4 ) (6) 



The field equation for a radially symmetric field is 

2(f) 3 - 2(f) 2 + 35(f)' 2 = 0, (7) 



d 2 o 2 do ., 2 



dr 2 r dr 
with the boundary conditions 

d(f) 

(f> — > (f) + as r — > oo, — = at r = 0. (8) 

ar 

Using the expression of we can calculate the Euclidean action of the bounce, which is given 
by 



Se(4>) = I d^-d^d^ + um 



4- jT rfr r 2 (i(^) 2 + [i(0 2 - I) 2 - - I) 2 ] + 5((f) 3 - 3 + )). (9) 



3 The bounce solution 

Before applying the systematic approach to solve Eq. (7), let us consider the thin wall 
approximation. This usually happens when 5 is small, and gives a solution nearly equals to 
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— 1 inside a sphere of radius R and nearly equal to +1 outside it. The region where <p differs 
significantly from these values is a kink of the from 

0(r) = tanh(r-i2), (10) 

which is the solution of the field equation 

d 2 (f> 2d<t> 3 2 



, • -2,r • 2,r 0. (11) 
dr A r dr 

One can easily show that the energy of the bubble (the action of the bounce solution) of 
radius R in the thin wall approximation (see for example [2]) is given by 

S E (R) = -^7iR 3 e + 4nR 2 S 1 , (12) 
where e — 2f and Si is the bubble-wall surface energy (surface tension), and is given by 

Si = jdriY^-Y + UM)}- (13) 

The critical bubble radius R, for which the energy is stationary, is written in terms of Si 
and e as 



whence it follows that 



16nSf _ 256tt 1 

5b - - "stf (15) 



Note that for in the thin wall approximation (i.e. small 5), we have exact analytical 
result for the radius as well as for the action. The radius is large and therefore the wall is 
indeed thin compared to the size of the bubble. 

For finite 8, the solution of the equation of motion Eq. (7) can't be written in a closed 
from. Following the approach in [5], the solution is constructed by means of an expansion 
in powers of 5. 

Since the bounce-solution is centered around the radius R, then it is more convenient to 
write the equation of motion in terms £ = r — R. Hence Eq. (7) becomes 

d 2 <p 2 d(p i 3 2 



v ^ + /^+ 2 * =° (16 > 

Depending on the thin wall approximation, we may write a Laurent series as an ansatz 
for the critical radius, 



4 



(I 

R = — + a + ai8 + a 2 5 2 + ... = V a^ . (17) 
The factor of the first derivative in the equation of motion becomes 

1 1 6 £ + ao x 2 , (g + Qo) 2 - Q i a -i ^3 , n / A 4N ,, fi N 

eTTTe + E^^ — ^r 5 + ^ 6 +0{6) - (18) 

Hence the expansion of the field equation in powers of 5 reads 

d 2 (f) 2 £ d(j) i 2(£ + a ) 2 d0 3 2 , 3 , 



+ 2 ^ + 2 r - 20 + 350 2 + 0(0 = 0. (19) 
at, 2 a_i a£ ar^ a£ 

Now, the solution of the above equation is obtained perturbatively up to the second order 
by means of the expansion 

= o + <50 1 + <5 2 2 + O(5 3 ). (20) 
For the zero order of S the equation of motion is 

-^ + 2^-20 o = O, (21) 

which has kink's solution 

MO = tanh£. (22) 

Note that the solution satisfies the boundary conditions given in Eq. (8), i.e., 0o(oo) — +1 
and <f>' (0) = 0. 

The equation of motion for the first order of 5 

-^2r - (6 sech 2 e - 4)0x = (— + 3) sech 2 £ - 3. (23) 
at, a_i 

The homogenous part of Eq. (23) has the following two independent solutions 

01,^(0 = -C lS ech 2 £ 

3 3 
0i,h 2 (O = C 2 (-£sech 2 £ + -tanh£ + sinh£cosh£). 

The particular solution of the equation which satisfies the boundary conditions for £ — > oo 
is 

0i,5 = C 3 + C 4 sech 2 £ 
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Substituting in Eq. (23), we get 

a-i = |, (24) 
which fixes the leading coefficient in R, Eq. (17), and 

<f>i,S = ~. (25) 
The general bounce-solution is given in the first order of 5 by 

faffi = -d sech 2 £ + C 2 & sech 2 £ + ^tanh£ + sinh£ cosh£) - -. (26) 

Z Z J_ 

Note that the constant term reflects the shift of the minimum given in Eq. (6). The term 
proportional to C2 diverges for large £, hence C2 must be zero. The first term corresponds to 
the derivative of the zero-order of the bounce-solution, and is related its translation degree 
of freedom. To see this, we expand the zero-order solution around 5 = in a Taylor series 

as 

tanh(£ - Ci<J) = tanh£ - d 5 sech 2 £ - C 2 5 2 sech 3 £ sinh£ + 0(5 3 ). (27) 

Note that C\ corresponds to the parameter a\ in the ansatz for the radius Eq. (17). The 
homogenous solution proportional to C\ is already taken into consideration by this ansatz, 
and should not consider in the the following calculations. Hence, we remain with 

MO = -| (28) 

Note that 4>i satisfies the boundary conditions given in Eq. (8). 
The equation of motion for the second order of 5 

d 2 (ho 9 3 9 

+ (6sech 2 £ - 2)0 2 = — £ sech 2 £ - - a sech 2 £ + - tanh£, (29) 



^ 2 v „ 8 ^ * 2 u * 8 

which has a general solution 

3 3 3 3 15 

2 (O = -Di sech 2 £ + D 2 (- £ sech 2 ^ + - tanh£ + sinh£ cosh^ + £(— + — cosh 2 £ - — sech 2 ^ 

— Log coshf (— £ sech 2 £ + — tanhf + — sinhf coshf) + — tanhf — — sinhf cosh£ 
& ^ V 32 S s 32 s 16 ; 64 s 16 

+ a^- + 4sinh 2 etanh 2 + ^ sech 2 £ T(£), (30) 



where we define 



T(£) = /Vtanh£'d£'. (31) 
Jo 



The general solution has a converging terms as well as diverging terms for large values of £. 
As the same argument for <f>i, D± — 0. The term proportional to = is symmetric in £, 
while the diverging homogenous solution of the bounce is antisymmetric, so we set ao = 
and we assume D 2 = jq — Log2 in order to get rid of the diverging terms. Therefore the 
converging solution is given by 

MZ) = -^(tanhe-^ + ^^cosK-sinhO 2 - ^'sech^+^sech^ + ^tanhe 
Log(l + e" 2 *) (A £ sec h 2 e + ^ tanh£ + ^ sinh£ cosh^ + ^ sech 2 £ T(£). (32) 

In order check whether 2 is an acceptable solution or not, we find its values at a large values 
of £, and its derivative at £ = —R. It can be easily shown that 

Hi oo) = ^, (33) 

which reflects the shift of the minimum given in Eq. (6), and 

0' 2 (£ = ~R) = 0. (34) 

Hence 02 satisfies the boundary conditions given in Eq. (8), which means that it is an 
acceptable solution. 

Whereas the first order solution corresponds to the shift of the minimum and the critical 
radius, the second order solution describes a true deformations of the bubble. The boundary 
condition at r = 0, i.e. £ = —R, is fulfilled order by order in 5. For example, the leading 
order solution yields 

4' (-R)=e- 2 "(4 + O(6)), (35) 

which vanishes to all orders in 5. Similar observations hold in higher orders. 
The bounce solution up to the second order of S is 

0(0 = O + 501+5 2 02 + O(5 3 ). 

= tan£ - 5- A + 5 2 ( - A £(tanh£ - 1) + A e ( cosh £ _ s[nh tf - A f sech 2 e + |_ ^ 

+ j| tanh£ - Log(l + ^ 2C )(t| £ sech 2 £ + A tan h£ + A sin h£ cosh^ 

+ A sec h 2 £r(0)+O(5 3 ). (36) 
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4 The action and the critical radius 



With the expression for we can calculate the action of the bounce, which is given by Eq. (9) 

as 



S E (<f>) = 47r^drr 2 (^ 2 +[i(0 2 -l) 2 -i(0 + -l) 2 ]+5(0 3 -^)) 



L2 

= S Q + SS 1 + 5 2 S 2 + 0(5 3 ). 



(37) 



The integrands are centered around the critical radius r = R. The integration range in £ 
can be extended to the whole real axis [5]. Hence 



/oo / 1 r 1 1 -i\ 



= 2tt(-R 2 



12 - 2tt 2 



)■ 



(38) 



Si = 4tt|^^(£ + j R) 2 (0[ ) / 1 + 20o0i^-1) + (^-1)) 



-00 
4 
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2 • (2 -)R). 



(39) 



/oo . u u . 

^df (£ + i?) 2 (0? + 2 0(,0' 2 + 6 2 2 - 2 2 + 40 o 2 (0 2 - 1) - - + -(1 - <&)) 

(40) 



' —00 
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= 27r(-^i? 2 + 1.01528). 



So, all together, the action is given up to the second order of 5 as 
S E = 2tt 



_ 12 2?f2 + 5 ( _ 1^3 + (2 _ Z^ )jR ) + S 2( _ l R 2 + L01528 ) 



(41) 



3 9 v 3 3 v - 4 

To find an expression for the critical radius R, we minimize the action with respect to 
R, i.e., 



dS 



E 



dR 



0. 



(42) 



Explicit calculation of the coefficients leads to two more terms in the Laurent series in R: 

(43) 



so that the bubble is now completely determined to the second order. Hence, the action 
follows as 



n r 128 1 .8 2tt 2 ^. r2 



(44) 
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For the thin wall approximation, we have obtained an explicit values of the radius and 
the bounce which are given in Eq. (14) and Eq. (15) respectively, these are exactly the values 
of the leading terms in Eq. (43) and Eq. (44). 

To check the validity of our analytical results, we have compared our analytical results 
with the thin- wall approximation for different values of 5. Figure 1 shows a plot of Se/Stwa, 
we can see that as far as 5 < 0.1, the ratio is of order 1, and for 5 > 0.1, the ratio is less 
than 1. The result is physically expected, since the maximum value of the action is in the 
case of thin-wall approximation. Similar observation also for the radius as shown in figure 
2. We have also solved the equation of motion (Eq. (7)) numerically. Figure 3 shows the 
ratio Euclidean action of the numerical results to the analytical results. We notice from the 
figure that as far as 5 < 0.1, there is a good agreement between them, and for 5 > 0.1, our 
analytical results starts deviates from the numerical ones. 

5 Conclusion 

By expanding all of the quantities in power of 5, we found the bounce solution, the critical 
radius and the action analytically beyond the thin wall approximation. We also showed that 
the leading terms correspond to the thin wall approximation. 

The comparison between our analytical results and the numerical ones showed a good 
correlation up to 5 = 0.1, and after that the analytical results starts deviates form the 
numerical ones. Also our analytical results are consistent with the TWA results up to 
5 = 0.1. 

To calculate the nucleation rate, we have to find the fluctuation determinant in powers 
of 5, this will be determined in future work. 
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Figure 1: Ratio of the analytical action to the action in the thin- wall approximation: 

Se/Setwa vs $ 
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Figure 2: Ratio of the analytical radius of the bubble to the radius in the thin-wall approx- 
imation: R/Rtwa vs 5 
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Figure 3: Ratio of the numerical action to the analytical action: Senum./ 'SEAnaiy. vs 5 
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